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We study nonlinear transport through a classical ballistic system accounting for the Coulomb 
interaction between electrons. The joint effect of the applied bias V and magnetic field H on 
the electron trajectories results in a component of the non-linear current I{V,H) which lacks the 
H — > —H symmetry: 51 = adV^H. At zero temperature the magnitude of ad is of the same order 
as that arising from the quantum interference mechanism. At higher temperatures the classical 
mechanism is expected to dominate due to its relatively weak temperature dependence. 
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By Onsager symmetry Q the electric current through 
a two terminal device in the linear regime is an even 
function of the external magnetic field, H . The nonlinear 
current does not necessarily satisfy this symmetry, as was 
recently observed in experiments IS 14 IsUsI ■ In fact 

it is violated already in the second order in the external 
bias V: It was shown in Refs. Hi that the nonlinear 
current acquires an odd in H component, 



is independent of the interaction parameter /hvp (here 
e is the electron charge and the vp is the Fermi velocity) . 
At small /hvp this is the leading contribution to ad- 
For an open two-dimensional ballistic dot the magni- 
tude of the classical contribution is 
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51 = aV'^H. (1) 

The violation of the H — > —H symmetry in the nonlinear 
current is associated with electron-electron interaction. 

Electron-electron interactions result in inelastic scat- 
tering of electrons. The standard phase space argument 
shows that the inelastic scattering cross-section varies as 
e^, with e being the electron energy, thus contributing 
to the current only at order and higher. Therefore 
in order to evaluate the coefficient a in Eq. o one may 
neglect inelastic processes and treat electrons as mov- 
ing independently in the presence of an effective poten- 
tial which depends on V and H. Within such a treat- 
ment HH the bias induces an additional electron density 
that is linear in the current and has an odd in H com- 
ponent. In the presence of interactions the additional 
density changes the scattering potential, which results in 
the nonlinear current O- 

In this approximation electron motion is phase coher- 
ent. Therefore a part of the nonlinear current is sensitive 
to the electron wave interference. The corresponding con- 
tribution to the coefficient a in Eq. (HJ was studied in 
Refs. 0, H IE llE It is a mesoscopic quantity whose in- 
dependence arises from the interference pattern sensitiv- 
ity to the magnetic flux threading electron trajectories. 

In addition to affecting electron interference the mag- 
netic field also bends electron trajectories. The corre- 
sponding contribution, ad , to the coefficient a is the sub- 
ject of the present work. In order to evaluate it electron 
motion can be treated as classical. 

We focus on the effect of the long range part of the 
Coulomb interaction and show that its contribution to ad 



where d is the dot size, and n, ep and m are the electron 
density, Fermi energy and mass, respectively. In Eq. © 
arises not from the electron-electron interaction but 
from the coupling of electrons to the magnetic field H 
and potential V, and from the definition of current I. 

At zero temperature the interference contribution is 0| 

4 

aq ~ where p is the density of states at the 

Eermi level, j3 is the the dimensionless interaction con- 
stant and Et is the Thouless energy. In the ballistic case, 
Et ~ hvp/d, and for /3 ~ I the coefficients ad of Eq. (|3) 
and aq are of the same order. Since aq decays with tem- 
perature on the scale T ~ Ep H ^^^*1 is insensitive 
to T for T < ep we expect the classical contribution to 
Eq. to dominate for T > Ep im 

The classical and the quantum interference contribu- 
tions can be distinguished by their respective scales of the 
magnetic field dependence. The linear i7-dependence o 
holds as long as the cyclotron radius significantly exceeds 
the system size d. This limitation yields the magnetic 
field scale H*^ ~ mvpcj{ed). The interference contri- 
bution Hi is sensitive to the flux threading a typical 
electron trajectory. Assuming that electron motion is not 
chaotic we estimate the area of such a trajectory to be 
d^. Equating the flux piercing a typical trajectory to the 
flux quantum <I>o = hcj e, we find the characteristic field 
for the interference contribution is much smaller than the 
classical one, Hq ~ ~ (Xp/d)H*p with Xp being 

the Eermi wavelength. 

Similarly to the quantum interference contribution H 
13 the magnitude and the sign of the classical contribution 
depend on the sample geometry. 
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While the existence of the nonlinear current Eq. o 
in a classical system and the estimate, Eq. ©, are quite 
general, below we concentrate on a specific setup consist- 
ing of a ballistic point contact in a two-dimensional elec- 
tron gas (2DEG) with an adjacent “reflector”, see Fig.^ 
The reflector creates the spatial asymmetry necessary for 
“rectification” current ~ V^. (The linear electron trans- 
port in systems of this type was studied experimentally 
in Ref. Il 2l l Electron motion in such a system is non- 
chaotic. This is essential for the validity of our results 
and the estimate The case of chaotic classical motion 
is beyond the scope of this work. 

In the classical description of electron transport the lo- 
cal momentum distribution function of electrons, n(p, r), 
plays a central role. The current I across the contact is 

where j(r) is the current density and the vector dS is 
normal to the contact line. A voltage bias applied to the 
contact results in a non-equilibrium, anisotropic electron 
distribution n(r,p) yielding a finite current I. 

In order to determine the nonlinear I-V characteristic, 
we need to find the steady-state nonequilibrium distri- 
bution function n(r,p). This can be done by solving 
the classical equations of motion for electrons moving in 
the presence of a self-consistent electrostatic potential 
</f(r). The initial conditions for the electron dynamics 
correspond to the equilibrium distributions deep inside 
the leads with electrochemical potentials differing by eV. 
Conservation of electron energy enables us to express the 
electron distribution at any point r in the form 



n(r,p) 



f{€p- fio + e(j){r) - eV), p e L{r), , , 

f{ep- fio + e(p{r)), p e i?(r). 



Here £p = / (2m) is the electron kinetic energy, /tq is the 

equilibrium chemical potential, /(e) = [exp(e/T) -I- 1]“^ 
is the Fermi function, and we assumed that the volt- 
age eV is applied to the left lead. The solution of the 
electron equations of motion is encoded in the shapes of 
the complementary momentum space domains, L{r) and 
i?(r). Electrons with momenta in these domains arrive 
to r from the left and right lead, respectively. 

The electric potential </>(r) is determined from the Pois- 
son equation. Its solution is greatly simplified if the 
screening radius is short compared with the geometri- 
cal characteristics of the setup. Under that condition, 
the charge density is unchanged by the applied bias. A 
further simplification is possible if the width w of the de- 
pletion layers confining the 2DEG is small compared to 
the geometrical features of the system. The electron den- 
sity increases from zero in the depletion region to its bulk 
value no over a distance of order w, see Ref. ES Thus at 
small w we may assume a constant electron density, 

/(W <"> 




FIG. 1: Point contact, thick lines 1 and 2, with an ad- 
jacent reflector, 3. The reflector line is i/ = L — x with 
{L — d) /2 < X < (T + d)/2. (a): Separatrix trajectories and 
the nonequilibrium electron distribution at R = 0. Vectors 
represent velocity directions on the separatrix trajecto- 
ries. Angles between and the x-axis are denoted by 
respectively. ( 6 ): Original (dashed line) and modihed (dash- 
dotted line) separatrix trajectories at R 7 ^ 0 , and a separatrix 
at R = 0 (solid line). 



inside the 2DEG. This is a reasonable approximation for 
devices fabricated using the local oxidation method Hi- 
The latter enables fabrication of structures with very 
small lateral depletion widths, w ~ 150 A ~ Af M- 

In order to find the current with the accuracy oc V^, 
we will solve the transport problem defined by Eqs. 

0 using the following iterative procedure. At the first 
step, we find the electron trajectories at />(r) = 0. This 
defines the zeroth iteration for the momentum domains, 
L^^\r) and The distribution function in this 

approximation depends on (jj{r) only explicitly, through 
the arguments of the equilibrium distribution functions 
in Eq. 0. Substituting n(p,r) into Eq. © we obtain 
an equation for (r) which is valid to the first order in 
V. After finding (r), we determine corrections to the 
electron trajectories and the corresponding corrections 
to the momentum domains. The corrected momentum 
domains L(r) and i?(r) and determine the distri- 

bution function via Eq. 0 . Substitution of the latter 
into Eq. 0 gives the current to the second order in V. 

For our device geometry each of the domains L{r) and 
R{r) consists of two separate sectors. The two sectors in 
L{r) correspond to trajectories arriving to r from the left 
lead either directly through the contact, or upon reflec- 
tion from the obstacle, see Fig.^ The boundaries of do- 
mains L{r) and i?(r) are defined by a set of four energy- 
dependent angles, 0^^^(r, Cp) between the x-axis and the 
velocity vector at point r for the separatrix trajectories. 
These trajectories pass through the edge points of the 
contact and through point r. The indices denote 
the trajectories that pass through the top(+)/bottom(“) 
edge point of the contact and arrive at r either directly 
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(d) or upon reflection off the obstacle (r). 

To evaluate the current I via Eq. 0 , we express the 
x-component of the current density, (r) , at the contact 
line, a; = 0, in terms of the angles e): 

jx(j) = 

Jxir) = 

Jxi^) = 
and treat I as the difference: 

I = Id-Ir, Id = j dSjiir ) , Ir = -J dSjliv). 

( 7 ) 

Here t'o = 2 ^p- is the electron density of states per spin 
projection, e_ = /Tq — e</)(r), e+ = /ig — e4>{r) + eE, and 
the integrals over dS go along the contact line. We set 
T = 0 for brevity. The generalization to finite temper- 
ature is straightforward; the scale for the temperature 
dependence of the current is set by the Fermi energy, /ig. 

Initiating the iterations, we find the angles 6^^^ = 

e) at V = 0. It is clear from Eq. 0 that to find 
the linear term in the I-V characteristic, it is sufficient 
to know (r, e) only at the Fermi energy. Equation 
© also indicates that the spatial dependence ^(r) does 
not affect the linear conductance M 

Next, using Eqs. 0 and © we obtain the first itera- 
tion for the potential 

^ (r) - (r) - (r) -k (r) 

V 27t 

( 8 ) 

Here the angles 9 are evaluated at the Fermi level, e = 
fio- This simplification is possible because the electron 
distributions in the domains L(r) and R{r) differ from 
each other only within a narrow energy strip of width eV 
around the Fermi energy. 

Finding the electric field from Eq. 0 we determine the 
correction to electron trajectories in the linear order in V. 
This defines the first iteration of the angles, 9^j^^\r,e). 
To find the current to the order we substitute 

e) = (r, e) + (r, /rg) (9) 

into Eqs. © and O- Note that we set e — > ^g in 9^^^ 
because 9^^^ is already proportional to V. 

In the absence of magnetic field and at zero bias the 
electron trajectories are straight lines. Therefore the an- 
gles 9^^^ do not depend on the energy and can be easily 
found from the geometric construction, see Fig. ^ 

= Imln[a; -I- iy T *d/2], (10a) 

^±,(o)(r) _ lniln[a: -I- iy — T(1 -I- i) ± d/2].(10b) 
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The first iteration 9^^'^ (r) is determined from the con- 
dition that the velocity direction at point r must change 
in the presence of the external force in such a way 

that the edge point of the contact still belongs to the sep- 
aratrix trajectory. A perturbative solution of Newton’s 
equations of motion gives, 

0(1) = r di f di'v (11) 

21qv^p m Jo Jo 

where /g is the length of the unperturbed trajectory, I 
is the coordinate along the trajectory, and Vj_Ji(i)(l') 
is the component of the electric potential gradient that 
is perpendicular to the trajectory. Substitution of the 
angles Eq. m found in the zeroth order into Eq. © and 
subsequent solution of Eq. m yield corrections 
which are proportional to V and depend on the geometry 
of the device. Substitution of the found 0^/°^ and 0^/^^ 

d/r d/r 

into Eq. yields the I-V characteristic with linear and 
cx terms: 



I = GshV 
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( 12 ) 



where Gsh = 2e^d.l {^\ f) is the Sharvin conductance of 
the point contact [l^. The rectification (oc F^) term and 
the corrections to the linear Sharvin conductance exist 
only due to the electron trajectories which are reflected 
by the obstacle back into the contact. The relative weight 
of such trajectories, contributing to Eq. ins), is oc d/L. 
In writing Eq. m we took the limit d L, in order 
to simplify the result. In this limit the rectification term 
arises solely from the energy dependence of the velocity 
m in the integrand of Eq. ©. 

Now we evaluate the influence of the magnetic field, 
assumed to be perpendicular to the plane of motion. In 
this case it is convenient to evaluate the “back-current” 
Ir in Eq. 0 in a different way. Equations (I6cll and 
express the “back-current” in terms of the directions of 
velocities of electrons that are already reflected from the 
obstacle and head towards the contact. We may, instead, 
evaluate the back-current by accounting for those trajec- 
tories of electrons that head towards the obstacle and will 
subsequently scatter back in the contact. We introduce 
the angles 9^{r) between the x-axis and the velocity di- 
rection at point r for the modified separatrix trajectories 
that start from r and, upon reflection from the barrier, 
arrive to the top/bottom edge point of the contact, see 
Fig. ^b). In terms of these angles, the back-current is 
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(13) 



Here the integration is performed along the contact line. 
The stationary-state current conservation law dictates 
/' = Ir- It is convenient for us to replace A in Eq. o 
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by {Ir +/')/2. Upon this substitution, we find: 



I 
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e/[r,ey 
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e+ (r,e) 
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(14) 

To arrive at this equation we used the fact that on the 
contact line 6*^(r) = ±7 t/ 2 + 0(iJ, V). This enabled us 
to set sin0^(r, e) = =pl in Efi- Jfibll . 

The modified separatrix trajectories coincide with the 
time-reversed original separatrix trajectories in the op- 
posite magnetic field and unchanged electric potential. 
Therefore the modified angles can be expressed in terms 
of the original ones as Of{r, H) = —H) — tt. Thus 

Eq. ini explicitly shows the H ^ —H symmetry of the 
linear current, in agreement with the Onsager relations. 

In the presence of a magnetic field, the trajectories be- 
come curved even at zero bias. We denote the deviation 
of the angle 9 from its value at = 0 by S9. For a 
trajectory originating from and arriving at without 
reflection off the barrier in a weak magnetic field we have 



S9^\rf-r,\/{2R,), (15) 



where Rc = mvc/{eH) is the cyclotron radius. The an- 
gles ( 2 :) acquire the correction 59'^^^\r) = \x + iy^ 
*d/2|/(2i?c). Though the expressions for the correspond- 
ing corrections to the reflected angles are more 

cumbersome they can be straightforwardly obtained from 
geometric considerations using Eq. Substituting 

these corrections into Eq. m we obtain the correction 
to the induced potential due to the magnetic field. It has 
an especially simple form on the contact line, x = 0, 



= o,y) 



u M' 

2uH = 

47ri?c [ ^ 2^2. 



(16) 



Although the second term here is caused by the reflected 
trajectories, the length L defining the barrier position 
drops out for L ^ d. Indeed, from Eq. 03 it follows that 
the corrections to the reflected angles are S9^ ~ L/Rc, 
whereas their difference entering the potential via Eq. (jS]) 
has an additional smallness of d/L. 

Next we evaluate the current across the contact using 
Eq. 11411 . It depends on magnetic field due to the H- 
dependence of the integration limits, e± , as described by 
Eq. itTHl) . and of the angles 9, 9 entering in the integrand. 
The finite- corrections to the angles 9 and 9 arise from 
two effects: i) the direct bending of electron trajectories 
by the Lorentz force, and ii) the iJ-dependence of the 
electric potential, Eq. itTCl) . For our geometry these ef- 
fects are small: The corrections to 9 and 9 due to the 
first effect are opposites of each other, as illustrated in 
Fig. G]b), and drop out from Eq. (tCTl . The correction 
due to the second effect arises only in the first iteration, 
and can be estimated as This turns out 

to be smaller than the correction arising from change of 



the energy integration limits e± in Eq. by a factor 
of d/L. Therefore, in order to find the coefficient ad of 
the V^H term in the nonlinear current to the leading or- 
der in d/L we may replace the argument of the sines in 
Eq. (inil by their H = 0 values. Doing so we obtain in 
the leading order in d/L 
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in agreement with the estimate Eq. @. 

While the “rigid-wall” confinement model applies to 
some setups 0 , the use of gates and scanning probes re- 
sults in a “soft” confining potential changing over length 
scales of the order of the size of the contact E3 This 
difference affects only the numerical coefficient in ad- 
In conclusion, we have shown that the nonlinear I-V 
characteristic of a classical ballistic system of interact- 
ing electrons lacks the H ^ — H symmetry. The classi- 
cal contribution to the odd in H current, Eq. , stems 
from bending of electron trajectories by a magnetic field, 
rather than from the flux sensitivity of electron wave 
interference pattern 0, 0- For a ballistic structure in 
which electron motion is not chaotic the magnitudes of 
the classical and interference contributions for weak fields 
H and low temperatures T are of the same order. The 
characteristic scales of H- and T- dependence for the 
classical contribution, however, exceed significantly the 
corresponding scales for the interference one. 
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